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RESUMEN 

Se presenta una serie de simulacioncs de 7V-cuerpos destinadas a estudiar la 
dinamica de grupos pequenos de galaxias. En particular, se comparan los resultados 
obtenidos con las propicdades dinamicas de los grupos de Hickson. Se consideran 
condiciones iniciales partiendo de una 'expansion maxima' y de equilibrio virial, 
y no se considera un halo oscuro comun primigenio. Se encuentra muy buena 
concordancia con las propiedades de grupos pequenos de galaxias, y las propiedades 
de los grupos compactos de Hickson son reproducidas adecuadamente por aquellos 
sistemas que se encuentran en un estado avanzado de un colapso gravitacional. 
No se encuentra un problema de 'sobre-fusion' de las galaxias. Para grupos que 
comienzan en virial, se encuentra que una fraccion importante (~ 40%) pueden 
durar por ~ 10 Gyr sin colapsarse completamente. Los resultados encontrados 
proporcionan otra solution alternativa al problema de sobre-fusion para los grupos 
de Hickson. Asimismo, se encuentra que la razon masa-luminosidad de los grupos 
de Hickson podrfa ser similar a la de cumulos de galaxias, sugiriendo que ambos 
tienen aproximadamcnte la misma fraccion de materia barionica a masa total. 

ABSTRACT 

A series of iV-body simulations aimed to study the dynamics of small groups of 
galaxies are presented. In particular, our results are compared with the dynamical 
properties of Hickson's compact groups (HCG). 'Maximum expansion' and virial 
initial conditions are tested, and no primordial common dark halo is considered. The 
properties of small galaxy groups are very well reproduced, and those of Hickson's 
groups are well reproduced by the most advanced stage of collapsing groups. We 
find no ovcrmcrging problem in our simulations. An important fraction of groups 
(~ 40%) initially in virial equilibrium can last for ~ 10 Gyr without complete 
merging. These results provide an alternative solution to the overmerging expected 
in Hickson's compact groups. Also, the mass-to-light ratio of HCG are probably 
similar to those found in clusters, suggesting that both kinds of systems have about 
the same fraction of barionic to total mass. 

Key Words: GALAXIES: INTERACTIONS - GALAXIES: KINEMAT- 
ICS AND DYNAMICS - METHODS: NUMERICAL 



1. INTRODUCTION 

Small groups are the most common galaxy as- 
sociations and contain about 50% of all galaxies in 
the universe (Huchra & Geller 1982; Nolthenius & 
White 1987). Important work has been done to com- 
pile catalogs with relevant kincmatical data. Thus 
Nolthenius & White (1987) and Nolthenius (1993) 
found that small groups from the CfA catalog (here- 
after NCfA) have the following median values: a one- 
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dimensional velocity dispersion of a m 116 km s , 
a mean harmonic radius of Rn = 480 kpc, a depro- 
jected median radius of i?s = 720 kpc and a crossing 
time of Hqt c = 0.44. 2 Hereafter a Hubble constant 
of Hq = 75 km s _1 Mpc _1 is assumed. Gourgoul- 
hon, Chamaraux & Fouquc (1992, hereafter GCF) 
found for nearby small groups (< 80 Mpc) a median 

2 Here, R s = 4R/n, R" 1 = 4£ i £ i< j Rr} / (7rN s (N g - 1)), 

and t c = 2iJfj/(\/3o"), where R, Rij and N s are the average 
projected separation, the projected separation between galaxy 
pairs and the number of galaxies, respectively. 
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o ~ 75 km s -1 , a virial radius of about 770 kpc 
and a crossing time-scale of ~ 4 Gyr, these quanti- 
ties corrected for the effect of the Hubble flow. More 
recently, Makarov & Karachentsev (2000, hereafter 
MK) derived from a sample of 839 groups, in particu- 
lar for those consisting of five galaxies, the following 
mean values, with quartiles: a = 57+^ km s -1 , a 
non-deprojected harmonic radius of 230tgg kpc and 
a dimensionless crossing time of about 0.1. 

A special and very important case is Hickson's 
compact groups (HCG) which are characterized for 
being composed of at least four galaxies with a me- 
dian projected separation of about three times the 
diameter of their luminous component and densities 
similar to the central region of rich clusters (Hick- 
son 1982; Hickson 1997). Hickson's initial catalog 
contained 100 groups with a median projected sep- 
aration of R = 52l 8 i kpc, a velocity dispersion of 
a = 200f3 7 km s~ x and a dimensionless crossing time 
of H Q t c = 0.016§;g 01 ; where t c = AR/^VSa), and 
subscripts and superscripts indicate the range of ob- 
served values. However, Sulentic (1997) pointed out 
that only 61 of them show according rcdshifts. 

It is noteworthy that HCG provide an ideal site 
to study the evolution of galaxies, via interactions 
and/or mergers, in a high density environment. One 
of the most striking challenges of the compact groups 
is how to reconcile their apparent longevity with 
their quite small crossing time of about 1% of the 
Hubble time. This short time-scale suggests a rapid 
evolution of these systems toward a complete merg- 
ing: the so-called overmerging problem (Hickson, 
Richstone & Turner 1977; Hickson 1982; White 1990; 
Hickson 1997). However, Hickson et al. (1984) and 
Rubin, Hunter & Ford (1991) have suggested that 
compact groups may have formed relatively recently. 

In general, two complementary approaches have 
been proposed to address this issue; see Athanas- 
soula (2000) for review on this topic. On the one 
hand, compact groups are assumed to be formed in 
the early universe and survive somehow until the 
present epoch (Ishizawa 1986; Diaferio, Geller & 
Ramella 1994; Governato, Tozzi & Cavaliere 1996). 
Diaferio et al. (1994) introduced a scenario where a 
continuous replenishment of the compact group oc- 
curs within a single collapsing rich loose group. Gov- 
ernato et al. (1996) proposed an alternative model 
that relies in the evolution of a compact group as 
an ongoing and frequent process through secondary 
infall of galaxies in a critical universe. However, 
in these scenarios merging is not suppressed so one 
would expect more compact groups and merger rem- 
nants than those observed today (Mamon 2000). 



On the other hand, there is a dynamical view 
where some kind of tuning of the initial conditions 
is imposed in order to guarantee the survival of the 
compact group. In this sense, Barnes (1985) was the 
first to investigate the effects of different kinds of ini- 
tial conditions on the evolution of compact groups. 
He found that (1) a compact group immersed in a 
massive common halo would, in general, delay its 
merging process if more mass is initially placed in 
this common halo. This merging delay results from 
the fact that the halo common mass is obtained from 
the galaxy haloes, thus reducing their masses and 
increasing the dynamical friction time-scale. This 
result was also corroborated by Bode, Cohen and 
Lugger (1993). And, (2), groups whose initial con- 
ditions begin in 'turnaround' or 'expansion' merge 
more rapidly than their virialized counterparts. 

More recently, Athanassoula et al. (1997) carried 
out a more systematic study of the parameter space 
for small groups (e.g. concentrations of the common 
and individual haloes, common-halo-mass to total- 
group-mass ratios, etc.). Their main findings are: 
(1) centrally concentrated groups merge faster for 
high common-halo to total-mass ratios and slower 
for low common-halo to total-mass ratios. (2) The 
overmerging issue of compact groups may not be a 
problem if appropriate initial conditions are chosen. 
To support this last point, they built up a virialized 
group with a large common-halo to total-group and 
an almost homogeneous central concentration that 
was able to survive for as long as <~ 20 Gyr. Ac- 
cording to these results the merger rates obtained 
from group simulations without including a common 
envelope, as the ones to be considered here, would 
be an upper limit. 

In this work the hypothesis of Barnes (1989) and 
White (1990) regarding that diffuse groups might be 
the progenitors of compact ones, but criticized by 
Diaferio et al. (1994) for lacking quantitative sup- 
port, is tested; similarly, the hypothesis advanced 
on observational grounds that HCG might be rel- 
atively young systems (Hickson et al. 1984, Ru- 
bin et al. 1991). The dynamical approach is used, 
and for this end a series of N-body simulations are 
performed without including a common halo. These 
simulations allow us to compare with the observed 
kinematical properties of small groups. The rest of 
this paper has been organized as follows: in section 
§2 the numerical methods to set up the galaxy mod- 
els and the initial condition for the group are de- 
scribed. In section §3 the merging histories of our 
simulations and the dynamical properties are given 
in section §4. In section §5 the mass estimation of our 
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groups is addressed. A constrained representation of 
the 'phase-space' is shown in section §6. A more re- 
alistic N-body simulation involving spiral galaxies is 
given in section §7, and section §8 contains a general 
discussion of the results. Finally, the main conclu- 
sions are summarized in section §9. 

2. DESCRIPTION OF THE SIMULATIONS 

This section contains a detailed description about 
the method used to set up our galaxy model and the 
initial conditions for our small galaxy groups. Also, 
the computational tools and the parameters used to 
evolve our simulations are described. 

2.1. Galaxy Model 

Individual galaxies are represented by a self- 
consistent Plummcr model. In building up the 
galaxy model no explicit difference between dark and 
luminous matter is done. However, to gain some feel- 
ing about the evolution of a luminous component, 
the initial 10% of most bounded particles are identi- 
fied as 'luminous'. 

By using the procedure described by Aarscth, 
Hcnon and Wielen (1974), positions and velocities 
of the particles are randomly generated from the fol- 
lowing mass and phase-space distributions /(£): 



M(r) = M g 



(r/Rof 



[1 + (r/i? ) 2 ]3/ 2 



and 



77r 3 G 5 M g 4 1 



where M g is the total galaxy mass, i? its scale- 
length, and £ the energy per unit mass, respectively. 

A system of model units such that G — 1 , M g = 1 
and Ro = 1 is chosen. Since at large radii the 
Plummcr density goes as p oc r~ 5 a cut-off ra- 
dius is introduced, for numerical purposes, at about 
r cut rj 10 i? which contains about 99% of the total 
galaxy mass. In this system of units, the dynam- 
ical time-scale, the half-mass and the cut-off radii 
are given by t d = ^8Rl/(GM g ) = VE, i? h « 1.3 
and r cu t = 10, respectively. In order to compare 
with observational data the values r cu t ~ 135 kpc 
and M g w 5.5 x 10 11 M© (Model B of Kuijken 
& Dubinski (1995)) are adopted. For these values 
the time and velocity units are 31.5 Myr and 419 



km 



In general, the simulations can be easily 



scaled through the following expressions for the ve- 
locity and time: v = 2.1 x 10 _3 ^M g /i? km s _1 and 
t = 4.7 x 10 5 ^/R^JW g Myr, where R and M g must 
be given in kpc and M©, respectively. With these 



definitions our initial 'luminous' core is contained 
within a radius of 0.74 model units, w 10 kpc. Fi- 
nally, each galaxy is represented by 3000 equal-mass 
particles. 

2.2. Group Initial Conditions 

In all cases each group consists of five equal- 
mass galaxies. Galaxies are assumed to be already 
formed and no secondary infall is considered, which 
is justifiable in a low density universe (e.g. Bah- 
call 1999, Hradecky et al. 2000). The positions and 
velocities for the center of mass of these galaxies 
are randomly generated from initial conditions start- 
ing at 'turnaround' and at virial equilibrium. In all 
simulations no primordial common dark halo is in- 
cluded. Although this assumption seems to contra- 
dict current cosmological models of structure forma- 
tion, some observations on clusters and groups sug- 
gest that most of the dark matter is associated with 
the individual dark haloes of galaxies (e.g. Bahcall, 
Lubin & Droman 1995, Puche & Carignan 1991). 
Hence, the merging activity observed in our simula- 
tions could be an upper limit. 

2.2.1. Collapsing groups 

For groups initially in 'turnaround' the positions 
for the centers of galaxies were obtained randomly 
from an homogeneous sphere of radius i? max de- 
fined by (Gunn & Gott 1972, ignoring any con- 
tribution from the cosmological constant) i? max = 
(8GMq£ 2 /7T 2 ) 1 / 3 , where M G represents the mass of 
the group enclosed within this radius at time t; mass 
conservation will be assumed. Compact groups are 
assumed to be at the verge of complete collapse at 
the present epoch (e.g. Hickson et al. 1984, Rubin et 
al. 1991), then their collapsing time-scale i c i ps turns 
to be the age of the universe to. Hence, i? ma x was 
evaluated at t /2 (Gott & Rees 1975). 

For fio = 1 it is found that to = 2i?g~ /3, and 
for n -» that t w H ( j\l + 0.5ft lnfi ) (Gott et 
al. 1974). 3 In particular for Q = 0.2 and H = 75 
km s _1 Mpc _1 to s=s 10 Gyr the 'turnaround' epoch 
was 5 Gyr ago (i.e. at redshift z w 0.6, Sandagc 
1961). In Table 1 the 'turnaround' radius of the 
group is listed for different values of Ho and Oo- A 
i? max = 700 kpc is taken as the fiducial 'turnaround' 
radius for collapsing groups. Notice that this value is 
consistent with a Ho = 50 km s _1 Mpc _1 and CIq = 1, 
and with the observational data of Sulentic (1987) 
who found that there is not a significant population 



3 This approximation differs by less than 1% from its true 
value for Q = 0.2 (Peacock 1999). 
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of normal galaxies within <~ 1 Mpc of the central 
region of HCG. 

Finally, the bulk velocity magnitudes for the 
galaxy centers were obtained randomly from a Gaus- 
sian distribution with a one-dimensional velocity dis- 
persion of cr = v /4GM G /(5i? max ) (Gott 1975) sub- 
ject to the constrain 2T/\W\ = 1/4 (Barnes 1985), 
where T and W arc the kinetic and potential energy 
of the group, respectively; these energies were com- 
puted assuming that galaxies are equal point-mass 
particles. Velocity vectors were oriented randomly. 

2.2.2. Virialized groups 

Theoretical arguments and A-body simulations 
(Peacock 1999) indicate that by a time w 3t c i ps /2 a 
system has reached an approximate equilibrium state 
with a virial radial scale of R v w i? max /2. Hence a 
radial scale for groups in near virial equilibrium at 
to is -R v = R(3to/2)/2, and in the particular case 
for our fiducial model a R v = 267 kpc is obtained. 
The positions for the galaxies are randomly chosen 
from a uniform spherical distribution with radius i? v , 
and their velocities from a Gaussian distribution re- 
stricted to satisfy the virial ratio 2T/\W\ = 1 with 
random orientation vectors. 



2.3. Computational Issues 

To evolve our galaxy and group models a serial 
TREECODE was used (Barnes & Hut 1986) with a 
tolerance parameter 8 = 0.9 and with the quadrupolc 
corrections to the potential included. A softening pa- 
rameter e — 0.07 w i?o/14 and a time-step At = 0.1 
were adopted. With these parameters, the simula- 
tions were evolved for about 321.4 time units (<~ 10 
Gyr) and the energy conservation, in all cases, was 
better than 0.7%. A total set of 30 simulations were 
performed for each kind of initial conditions consid- 
ered for the groups. 

To follow up the orbital evolution of each galaxy 
inside the group its center was identified with its 
most bounded particle, which was set by 'hand' to 
have 1% of the galaxy mass (Aguilar & White 1985). 
In order to check the stability of the galaxy model 
due to this change its evolution was followed, in iso- 
lation, for <~ 200 dynamical time-scales, no signifi- 
cant change in its inner structure and in the total 
energy conservation was found. From here on, all 
kinematical properties of the group simulations (e.g. 
mean separation, velocity dispersion, etc.) will be 
obtained from these galaxy centers. 



TABLE 1 
TURNAROUND RADIUS (MPC) 



Ho 


n = 0.2 


n = 1 


(km s^Mpc" 1 ) 






50 


0.88 


0.76 


75 


0.67 


0.58 


100 


0.56 


0.48 



3. MERGING HISTORIES 
In this section the merging histories for all our 
simulations of small groups are computed. These 
histories provide some insight about the merging 
rate and the ovcrmcrging problem of compact groups 
posed by their small crossing time. 

In Figure 1 the evolution of three particular sim- 
ulations with cold initial conditions is shown. No- 
tice that group g27c does not present any merger 
during the simulation (~ 10 Gyr) while group g02c 
completely merged. At present epoch (i.e. ~ 5 Gyr 
from 'turnaround'), group g22c resembles a chain- 
like compact group which is better appreciated at the 
bottom-right of this figure where only 'luminous' par- 
ticles have been displayed. Notice also that groups 
<?27c and g02c may be classified as triplets. 

In order to characterize this merging activity the 
following pair-wise merging criteria (Aarseth & Fall 
1980) are employed: 

V ij = l V j - V *l < Vrms/2, Rij = \rj - rt\ < Rq/2 ; 

where Vij and are the relative velocity and sepa- 
ration between a pair of galaxies i and j, respectively, 
and V rms and Rq are their characteristic internal ve- 
locity and radius. 

In Figures 2 and 3 the time evolution of the 
merging activity in our groups is shown. The 
dynamical time for the virialised groups is td — 
y/8R^/(GM(j) i=a 111 corresponding to about 3.6 
Gyr. It is observed in Fig. 2 that no group has 
merged completely at the present epoch, to, although 
some multiple mergers do exist. About 15 groups 
(50%) do not present any merger at all, 9 (30%) 
show a merger of two galaxies and 6 (20%) triple 
or quadruple mergers. Notice that at the end of the 
simulations just one group does not show any merger 
and 15 (50%) have merged completely. 

In the case of virialized groups, 5 (~ 17%) do 
not present any kind of merging within a dynamical 
time-scale, 11 (ss 37%) show a merger of two galax- 
ies, 10 (w 33%) have triple or quadruple mergers, 
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Fig. 1. This figure shows the XY-projection in kpc for the collapsing groups (/27c (upper-left), <?02c (upper-right) and 
i?22c (bottom-left). The numbers in each frame indicate the time elapsed since turnaround in Gyr. Lines refer to the 
trajectories followed by each galaxy inside its corresponding group. Also, small boxes of 300 kpc wide are shown at time 
4.79 Gyr for each group. These boxes have been amplified and show just 'luminous' particles at the bottom-right of this 
figure which corresponds to the XY, XZ and YZ-projections. 
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Fig. 2. Number of individual galaxies versus time 
(merging history) for simulations with cold initial con- 
ditions. The free-fall time-scale (Mamon 1990) tg = 
it^J -R? m ax/(32GMG) and the present epoch to (~ 5 Gyr 
from 'turnaround') are indicated by arrows. The merging 
tracks have been displaced vertically and in time for bet- 
ter appreciation. The displacement in time corresponds 
to 0, 3, 6, 9 and 12 time units. A scale corresponding to 
1 Gyr is indicated. 




100 
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Fig. 3. As in Fig. 2, but for initially virialised groups. In 
this case, the arrow corresponds to the dynamical time 
of the group. 
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and only 4 (i=s 13%) have completely merged (see 
figure 3). Notice that about 12 (40%) have not com- 
pletely merged after ~ 10 Gyr of evolution. This 
clearly indicates that some compact groups can sur- 
vive complete merging for about a Hubble time. 

4. DYNAMICAL PROPERTIES 

The dynamical properties of the groups along 
three orthogonal 'lines-of-sight' at different times 
are determined. The median values for the one- 
dimensional velocity dispersion, the mean harmonic 
and deprojected radii, the VME and MME masses 
(see §5), and the dimensionless crossing time, Hqt c , 
and the dimensionless crossing time as in Hickson 
ct al. (1992), t c , are summarized in Table 2. The 
times indicated in the table correspond to 0.5, 1, 2, 

10 Gyr of evolution from their respective initial 
conditions. 

4.1. Collapsing Groups 

In Figure 4 the dynamical parameters as a func- 
tion of time for our collapsing groups g27c, g02c and 
gTlc are shown, respectively. For comparison pur- 
poses the medians of MK and HCG are indicated by 
arrows. This figure shows the general trend in a col- 
lapse; a decrease in size followed by an increase in 
velocity dispersion and a reduction in the crossing 
time. The size of a group depends on the method 
adopted; the harmonic radius is always lower than 
the mean separation radius and, hence, r c < t c . An 
agreement between both definitions of crossing times 
occurs when the systems are still rather diffuse. This 
is due to the fact that the harmonic radius i?H is 
more sensitive to smaller separations than i?g- 

The kinematical properties of g21c and g02c at 
the present epoch are in good agreement with the 
values found in NCfA, GCF, and MK groups. In gen- 
eral, our simulations reproduce very well the kine- 
matical parameters of 'normal' small groups of galax- 
ies (see Table 2). 

Group g27c. This group has properties at the 
present epoch that resemble those of a 'normal' dif- 
fuse group: i?H ~ Rs ~ 700 kpc, a w 70 km s _1 , 
and Hqt c w H^t c w 0.5. These values are large with 
respect the observed ones in HCG, but in better 
agreement with the medians of NCfA, GCF, and MK 
groups. 

Group g02c. This system has a similarity with 
Hickson's groups; it shows a a w 200 km s" 1 and 
Hot c w 0.09 when is observed along the Y-axis (see 
Fig. 4). However, its mean separation is quite large 
Rs ~ 400 kpc. The average values, over the three 
lines of sight, obtained at the present epoch are: 



1 

0.8 
0.6 
0.4 
0.2 




Gyr 



-1 
Log fj 





-10 12 
Log R H 



1 rrr 

0.8 - 

0.6 - 

0.4 - 

0.2 - 
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Fig. 5. Collapsing groups: Histograms for the velocity 
dispersion, mean harmonic radius, dimensionless cross- 
ing time (Hot c ), and mean separation at different times. 
The different symbols denote the median for HCG (*), 
NCfA (|), and MK groups (j). Hickson Compact Groups 
Associations' values (Ramella et al. 1994) are also indi- 
cated (J.) for completeness. Different line styles corre- 
spond to different times of evolution: 0.5 Gyr (dotted), 
3.0 Gyr (short-dashed), and 5 Gyr (solid). 



R H « 10 kpc, R s w 250 kpc, er « 90 km s" 1 , 
H t c w 0.01, and H t c w 0.2. 

Group g22c. This group has dynamical proper- 
ties, at the present epoch, that resemble very closely 
the medians of Hickson's groups: i?H ~ 40 kpc, 
R s « 110 kpc, a w 120 km s -1 , H t c « 0.1, 
and Hot c w 0.05. However, a Rs ~ 50 kpc, 
a w 180 km s _1 and Hot c ~ 0.01 are found along 
the Z-axis. We should mention that the selection of 
g22c was based solely on its similarity with HCG at 
the present time when seen in the XY-plane, which 
is confirmed with the more realistic models of spiral 
galaxies in section §6. In this group, the three closest 
galaxies merge in ~ 1 Gyr from the present epoch, 
while the other two form a wide 'binary' at the end 
of the simulation (10 Gyr). This result suggests that 
although some HCG have small crossing times they 
will not merge completely over several Gyr from the 
present epoch. 

Histograms for the dynamical properties of col- 
lapsing groups at times t « 0.5, 3, 5 Gyr, normalized 
to the initial number of 'observed' groups, are shown 
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TABLE 2 

DYNAMICAL PROPERTIES OF GROUPS IN MODEL UNITS 









Collapsing 














Virialized 








t 


G 


Rh 


H T C 


Rs 


H t c 


M v 


M mcd 


a 


Rh 


H T C 


Rs 


H t c 


M v 


M med 


15.0 


0.071 


45.76 


1.91 


51.27 


1.09 


1.28 


0.65 


0.221 


15.14 


0.197 


19.17 


0.126 


3.41 


1.79 


36.0 


0.077 


43.25 


1.57 


49.74 


0.99 


1.50 


0.67 


0.223 


10.79 


0.126 


17.18 


0.111 


2.71 


2.06 


63.0 


0.090 


36.09 


1.01 


46.78 


0.72 


1.46 


0.79 


0.189 


2.72 


0.043 


18.65 


0.134 


0.99 


1.49 


96.0 


0.112 


25.65 


0.43 


41.19 


0.49 


1.33 


0.62 


0.175 


1.44 


0.033 


20.10 


0.179 


0.27 


0.66 


126.0 


0.132 


8.41 


0.11 


34.04 


0.33 


0.88 


0.92 


0.141 


0.55 


0.015 


21.23 


0.231 


0.06 


1.22 


156.0 


0.149 


3.70 


0.06 


27.86 


0.28 


0.76 


1.13 


0.137 


0.45 


0.012 


20.31 


0.212 


0.06 


0.60 


186.0 


0.153 


0.87 


0.02 


22.31 


0.21 


0.16 


1.74 


0.111 


0.24 


0.008 


21.47 


0.195 


0.02 


0.08 


219.0 


0.158 


0.53 


0.02 


17.21 


0.19 


0.11 


1.10 


0.111 


0.21 


0.008 


21.10 


0.241 


0.01 


0.04 


252.0 


0.152 


0.28 


0.01 


15.07 


0.17 


0.04 


0.95 


0.082 


0.15 


0.007 


20.66 


0.239 


0.01 


0.02 


282.0 


0.131 


0.23 


0.01 


15.59 


0.20 


0.03 


0.57 


0.077 


0.15 


0.007 


18.39 


0.242 


0.01 


0.01 


315.0 


0.100 


0.19 


0.01 


17.20 


0.21 


0.02 


0.14 


0.083 


0.11 


0.005 


15.41 


0.197 


0.01 


0.00 



t time evolved from the initial conditions. 

<r one- dimensional velocity dispersion. 

Rh harmonic radius. 

Rs mean separation radius. 

M v virial mass. 

M med median mass. 

H a t c crossing-time scale where t = R s /(VScr). 
H r c crossing-time scale where t c = 2_R h /(v / 3<t). 



in Figure 5. The median values obtained by differ- 
ent authors are also indicated. The crossing time 
as defined by Hickson et al. (1992) is shown in the 
left-panel of Figure 6. 

This ensemble of simulations shows a general 
trend to lower their size and increase their veloc- 
ity dispersion, until the point where 'binaries' or 
'triplets' start to form. At this stage the velocity 
dispersion of the group begins to decrease due to the 
transfer of orbital energy to the internal degrees of 
the individual galaxies. Also, the harmonic radius 
distribution shows a strong increase in its lower end 
while the velocity distribution increases its higher 
end (see Fig. 5). Notice that at the present epoch 
(i w 5 Gyr) there is a significant fraction of groups 
with a > 100 km s _1 . 

It is important to point out that the median val- 
ues observed in catalogs of 'normal' small groups 
are well reproduced by our collapsing groups at the 
present epoch (e.g. a and Rs in Figure 5). About 
30% of the collapsing groups have a a value close 
to the median of HCG and w 50% have a similar 
Rs- Furthermore, the majority of the crossing times 
fall within the values found in HCG (Figure 6) and 
w 50% would be classified to be virialized at the 
present epoch accordingly with the criterion of Gott 



& Turner (1977). 

Hickson's catalog (Hickson 1982, Hickson et al. 
1992) has 16 groups with a < 100 km s _1 and just 
five of them show a a < 5 km s _1 . Four of these last 
groups (HCG38, HCG47, HCG49 and HCG88) have 
a dimcnsionless crossing time of H n t c > 1 suggesting 
they are not close to virial equilibrium. However, our 
results indicates that such a assumption may not be 
correct. 

Since our models roughly agree with those HCG 
exhibiting a low velocity dispersion, they suggest 
a probable explanation for that fact (Rubin et al. 
1991, Mamon 2000). 

4.2. Virialized Groups 

Figure 7 contains histograms for the dynamical 
quantities for virialized initial conditions at times t ~ 
0.5,3,6 Gyr. According with Table 2, these groups 
present a median of a w 90 km s _1 , Rh ~ 200 kpc, 
R s « 250 kpc, H t c w 0.20 and H t c w 0.13 at their 
initial stages. After 5 Gyr of evolution the following 
medians are obtained: a ~ 58 km s _1 , Rn w 6 kpc, 
R s w 274 kpc, H t c w 0.01 and H t c w 0.21 which 
are far from the values found in Hickson's groups. 

In general, the cr-distribution shows small 
changes even after ~ 10 Gyr of evolution; however, 
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Fig. 4. Kinematical parameters as a function of time for collapsing groups g27, g02c, and g22c. Astronomical units 
are used. The left-arrow (<— ) indicates the median of the 5-galaxy MK groups, while the right-arrow (— >), HCG. The 
different lines indicate the line-of-sight used to compute the quantities: along X-axis (solid), along Y-axis (dotted), and 
along Z-axis (broken). For the dimensionless crossing times, thicker lines for Hqt c than for Hot c have been used. 



the i?u-distribution presents significant changes 
which are also observed in the evolution of the Hqt c - 
distribution. The tr-distribution tends to populate 
the regions of lower velocity dispersion and to in- 
crease somewhat the mean separation. Notice that 
about 50% of groups have at t = crossing times 
similar to those of normal small groups (NCfA, GCF, 
and MK) and ~ 5% close to the median of Hickson's 
groups. 



dicate that a large fraction (~ 40%) of our initially 
virialized groups have not merged during 10 Gyr and, 
hence, the overmerging (Hickson 1997) may not be 
a problem. However, it is important to point out 
that in our simulations a primordial common dark 
halo has not been included, which may change our 
conclusions. 



The results of this and the previous sections in- 
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Fig. 6. Evolution of the distribution of the crossing-time t c = Rs/(V3a). Lines styles are the same as in Fig. 5. The 
interval of crossing-times for HCG and that assumed to be satisfied by virialized groups, according the criterion of Gott 
& Turner (1977), are indicated. The left panel shows results for collapsing groups, while the right one for virialized. 
Notice that both kinds of initial conditions lead to ranges in crossing-time within the observed values. 
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Fig. 7. Same as in Fig. 5, but for small groups with virial 
initial conditions. 



5. MASS ESTIMATES 
The evolution of the virial mass (VME) and me- 
dian mass (MME) estimators is computed. These 
estimators are defined, respectively, as (Heisler, 
Tremaine & Bahcall 1985, Aceves & Perea 1999 and 
references therein): 



2G E^ l/fl 



M med = ^MED.y [(Vi - Vj) 2 Rij\ , 

where N s is the number of galaxies, Vi the velocity 
along the line-of-sight with respect to the centroid of 
velocities, Rij the projected separation on the sky, 
and the constant in M me d is determined from numer- 
ical experiments. 

In Figure 8 the above quantities are plotted as 
function of time for our collapsing groups <?27c, <?02c 
and g22c. This figure reveals that both mass estima- 
tors underestimate the mass of the group at t = 
as expected (remember that these numerical groups 
started with a virial coefficient far away from equi- 
librium 2T/\W\ = 1/4). As the group evolves and 
reaches some state of virial equilibrium it is observed 
that these mass estimators give a reliable mass for 
the groups as long as no mergers occur between their 
galaxies members. Notice that only bulk motions 
are considered in their definitions ignoring the self- 
gravitating nature of the individual galaxies. This 



last point is important since self-gravitating galax- 
ies have the ability to absorb energy as the whole 
group evolves. In general, the VME method is more 
sensitive to the dynamical state of the system than 
its mass, while the MME is dominated by the large 
median separations and low velocity dispersions. 

The results from all the simulations carried out 
are summarized in the histograms of Figure 9. The 
true mass of the groups is indicated by the vertical 
line. For collapsing groups (left panels) the VME es- 
timator (left- upper panel) gives a median virial mass 
of My = 1.28 at t = 15 which is consistent with the 
expected value of My = 1.25, 4 while the MME (left- 
lower panel) underestimates the mass by about an 
order of magnitude. However, once the groups are 
allowed to evolve to the present epoch (~ 5 Gyr) the 
mass is underestimated by an order of magnitude for 
the VME and by a factor of about 5 for the MME 
method. For the case of virialized groups we found 
that at t = both mass estimators give a similar 
median mass of about My w 3.5. Once the separa- 
tion between galaxies decreases, the MME provides a 
better estimate of the mass. From these results, the 
MME is considered more appropriate to estimate the 
mass of systems that resemble HCG. 

Hickson's compact groups have a mass-to-light 
ratio of T ss 50/iTq. If HCG are physically well de- 
fined groups, their content of matter is probably > 5 
times more than is inferred now. This would im- 
ply that a closer estimate to the true mass-to-light 
ratio of HCG would be T ~ 250/iY©. It is inter- 
esting to note that this value falls within the values 
found in clusters 200 - 300/iT Q (Bahcall 1999) sug- 
gesting that both types of systems would have ap- 
proximately the same fraction of baryonic to dark 
matter (see White et al. 1993 for a discussion of the 
cosmological implications) . 

6. R s ~a DIAGRAM 

A restricted 'phase-space' consisting of a i?s _ c 
diagram is constructed in order to try to discern 
some evolutionary track in our simulated groups. In 
Figure 10 the collapsing (left-panel) and virialized 
(center-panel) groups are shown at times t = 0.5 
and t — 5 Gyr together with the corresponding dia- 
gram for HCG where the ones considered 'bona fide' 
by Sulentic (1997) are indicated. 

This diagram clearly shows that initially cold 
groups reduce their mean separation and increase 



4 According to Bahcall & Tremaine (1981) and Geller 
(1984) a lower fractional limit on the virial mass is <5 V ^ 
v / 2mA7(7r 2 v / A f ), leading to M v ss 1.28 ± 0.12 for N = 90 
'observed' groups. 
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Fig. 8. Time behavior of the two mass estimators considered in §5 for our collapsing groups (/27c, (/02c and (/22c. The 
upper panels show the virial mass while the bottom the median mass. The horizontal broken line corresponds to the 
total mass of the groups; iV-body units are used for the mass while physical units for the time. The different types of 
line denote different lines-of-sight as in Figure 4. 
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Fig. 9. Histograms of the behavior of the VME (upper panels) and MME (lower panels), both for cold (left) and 
virial (right) small groups. The bar (|) indicates the true mass of the groups in iV-body units. Dotted lines show the 
distribution at t — 0.5 Gyr, short-dashed lines at 3 Gyr and solid at 5 Gyr. 
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their velocity dispersion as they evolve, while for 
virializcd groups their velocity dispersions decrease 
somewhat and their mean separation tends to in- 
crease. Furthermore, about 8 collapsing groups have 
R s < 100 kpc and a <= (60,200) km s" 1 at the 
present epoch (~ 5 Gyr) which roughly correspond 
to the values shown by 'normal' small groups, and 
are within the range of values found for HCG. Also, 
2 groups have R$ £ 50 kpc and a rts 200 km s _1 
which are very close to the median values of HCG. 

This Rs-p diagram suggests a resemblance of 
HCG with our the systems that once had a dif- 
fuse configuration and are at the verge of collapse 
at the present time, while some of our virialized 
groups show properties more similar to the 'normal' 
small groups. These and previous results suggest 
that HCG are the more advanced systems in a pro- 
cess of gravitational collapse among small groups. 
Thus, compact groups may be rather young config- 
urations (Hickson et al. 1984, Barnes 1989, White 
1990, Rubin ct al. 1991). 

7. A MORE REALISTIC SIMULATION 

In the simulations reported above no clear dis- 
tinction between luminous and dark matter was 
made, and hence no precise statement could be es- 
tablished whether, for example, the group g22c could 
actually be considered a compact group according to 
Hickson's criteria. 

To address this particular issue a more realistic 
galaxy model was set up by replacing the plummer 
models by spiral galaxies using model B of Kuijken 
& Dubinski (1995); see §2.1 for some numerical val- 
ues of this model (we refer the reader to Kuijken 
& Dubinski's paper for a detailed description of the 
method to build up the disk galaxy model). It should 
be mention that a detailed study of groups involving 
disk galaxies is out the scope of the present work. 
Table 3 summarizes the parameters that define our 
galaxy model and Figure 1 1 shows its corresponding 
rotation curve. 

For the galaxy centers and bulk motions inside 
the group the values of our previous collapsing model 
g22c are used. For illustrative purposes, arbitrary 
orientations for the galaxy angular momentum vec- 
tor are chosen. Each numerical galaxy consists of 
65536 particles, 16384 for the disk, 4096 for the bulge 
and 45056 for the halo; i.e., a total of 327680 were 
used in this simulation. 

To evolve the system a parallel treecode de- 
veloped by Dubinski (1996) with a tolerance param- 
eter of 0.8 was used. Forces were computed with the 
quadrupole terms included and with a fixed time- 
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Fig. 11. Rotation curve ftf k t>he spiral galaxy model. 



step corresponding to 1 Myr. With these parameters 
the energy conservation was better than 0.2%. The 
simulation was ran in a Beowulf class cluster consist- 
ing of 32 Pentium processors of 450 MHz (Velazquez 
& Aguilar 2002). It took a wall time of 16 seconds 
per time-step. 

In Figure 12 (left) the XY-projection of the lu- 
minous component of the group at time 4.8 Gyr, 
the same as figure 1 of model g22c, is shown. Al- 
though differences in the evolutionary trend of our 
more realistic model are appreciated with respect to 
model <?22c, it is rather easy to conclude that it sat- 
isfies Hickson's criteria for compact groups. Hick- 
son's criterion on the difference of magnitudes among 
galaxies is satisfied immediately since our galaxies 
are identical. Therefore, a group is obtained that 
is compact at the present epoch and will not merge 
completely for another ~ 5 Gyr. 

In Figure 12 (right) only the dark particles of the 
haloes are shown. Notice that a common dark halo 
has already formed for three of the galaxies that will 
merge in the next few Gyr. This result is consis- 
tent with the suggestion by Rubin ct al. (1991) that 
compact groups reside in a common dark halo at the 
present epoch. 

8. GENERAL DISCUSSION 

It has been found that w 10% of initially dif- 
fuse groups reach the present epoch with dynamical 
properties consistent with the median of Hickson's 
compact groups. Hickson's groups appear as those 
systems that have manage to evolve more rapidly 
toward a compact configuration, due to particular 
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Fig. 10. Rs-cr diag^alfflfor collapsing (left), virialized (center), and Hickson's groups (right). Hickson's compact groups 
appear to resemble systems that have suffered a gravitational collapse out of a more diffuse system. 



Fig. 12. (left) XY-projection at time 4.8 Gyr of our more realistic model. Luminous particles are shown in the left-panel 
and the dark component is indicated in the right-panel. Notice that a common dark envelope for the three galaxies close 
to merging is already present, while the other two in the upper-left will not merge for > 5 Gyr from the present epoch. 
Compare with our previous model p22c of Figure 1. 
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TABLE 3 
GALAXY PARAMETERS 



Model 


M D 


Rd 


Rt 


ZD 


SRout 




OB 


PS 


*o 


O"0 


q 


c 


Ra 


MW-B 


0.87 


1.0 


4.0 


0.15 


0.4 


-2.9 


0.71 


14.5 


-5.2 


0.96 


1.0 


0.1 


0.8 



M D , _R d , i? t , z D and <5_R ollt correspond to the mass, radial scale-length, cut-off radius, vertical scale- length of the disk, and 
disk truncation width respectively. * c , a B and p B refer to the cutoff potential, velocity dispersion, and central density 
of the bulge, respectively. Finally, * , cr , <?, C and R a indicate the central potential, velocity dispersion, potential 
flattening, concentration, and characteristic radius of the halo, respectively. 



initial conditions in the density perturbation that 
led to the formation of small groups. The forma- 
tion of a compact groups is, according to our re- 
sults, consistent with a hierarchical clustering sce- 
nario for the formation of structures in the uni- 
verse. Hence supporting the hypothesis advanced by 
Barnes (1989) and White (1990), and the one sug- 
gested on observational grounds (Hickson et al. 1984, 
Rubin et al. 1991). 

The present model leads to ~ 50% of originally 
diffuse groups arriving at the present epoch without 
any merging. Also, it is found that ~ 10% (5 sim- 
ulations) of these groups will not have any merger 
within the next <~ 5 Gyr. In this sense, our model 
provides an alternative explanation to the proper- 
ties of compact groups and their existence today, and 
complements several previous studies (e.g., Athanas- 
soula et al. 1997). 

In general, initial conditions from 'maximum ex- 
pansion' reproduce very well the properties of 'nor- 
mal' small groups. Our simple model avoids the re- 
quirement for an initial overdensity and/or signif- 
icant secondary-infall (Diaferio et al. 1994, Gover- 
nato et al. 1996), and the need for especial initial 
conditions to avoid the complete merging of a group 
(e.g., Athanassoula et al. 1997). 

The present model does not show the merging 
instability, but many questions remain. One of the 
most obvious caveats of our simulations is not taking 
into account a primordial dark halo for virializcd ini- 
tial conditions. However, the merging activity found 
may still be considered a lower limit if a concentrated 
dark halo were to be included. Surely more realistic 
simulations of small groups, considering a spectrum 
of masses and a clear distinction between luminous 
and dark matter, as well as initial conditions taken 
from a cosmological simulation, are necessary to ad- 
dress these issues. We plan to study some of these 
matters in the future. 



9. CONCLUSIONS 
The main results of this work are as follows: 

1. Groups starting from 'maximum expansion' 
with a diffuse configuration have, at the present 
epoch, dynamical properties very similar to the 
medians obtained from several catalogs of small 
groups. In particular, about 10% of them show 
values close to the observed medians of Hick- 
son's groups. 

2. The suggestion that Hickson's compacts groups 
originated from diffuse systems, and are rela- 
tively young systems, finds quantitative support 
in our numerical experiments. 

3. It is found that overmerging is not an impor- 
tant problem in our simulations, either with viri- 
alized or cold initial conditions. For the case 
of virialized groups about 40% survive merg- 
ing during « 10 Gyr. This also indicates that 
the existence of a large, massive common dark 
halo does not appear to be a strong requirement 
to explain the present day properties of small 
groups. 

4. The median mass estimator (MME) appears to 
be a better estimator of mass than the virial 
mass estimator (VME) for groups at the present 
epoch. It is estimated that the mass-to-light 
ratio of HCG is probably T ~ 250T Q suggesting 
that clusters and compact groups have about 
the same fraction of baryonic to dark matter. 
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